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the properties of latex particles forming a pressure sensitive adhesive. Each latex particle 
is treated as a single extensible particle. Transient forces are essential to correctly predict 
qualitative features. 
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Abstract 

For optimal application, pressure-sensitive adhesives must have rheological properties in between 
those of a viscoplastic solid and those of a viscoelastic liquid. Such adhesives can be produced by 
emulsion polymerisation, resulting in latex particles which are dispersed in water and contain long- 
chain acrylic polymers. When the emulsion is dried, the latex particles coalesce and an adhesive 
film is formed. The rheological properties of the dried samples are believed to be dominated 
by the interface regions between the original latex particles, but the relation between rheology 
and latex particle properties is poorly understood. In this paper we show that it is possible to 
describe the bulk rheology of a pressure-sensitive adhesive by means of a mesoscale simulation 
model. To reach experimental time and length scales, each latex particle is represented by just one 
simulated particle. The model is subjected to oscillatory shear flow and extensional flow. Simple 
order of magnitude estimates of the model parameters already lead to semi-quantitative agreement 
with experimental results. We show that inclusion of transient forces in the model, i.e. forces 
with memory of previous configurations, is essential to correctly predict the linear and nonlinear 
properties. 



Corresponding author, e-mail: j . t . paddingOgmail . com 



2 



I. INTRODUCTION 



Pressure sensitive adhesives (PSAs) are tacky viscoelastic or viscoplastic materials that 
adhere to a substrate upon the apphcation of hght pressure [l| . A well known example of a 
PSA is the sticky layer on Scotch@ tape. Many modern PSAs are soft nanostructured mate- 
rials formed by drying of an aqueous dispersion of polymeric particles, also known as a latex. 
Such polymeric adhesives offer many advantages over rubber-based or ceramic alternatives, 
like for example a low plateau modulus, superior stability, great resistance to oxidation, 
and avoidance of harmful or volatile organic solvents. Unfortunately, the design of new 
PSAs is hampered by a lack of fundamental understanding of the link between microscopic 
parameters, such as polymer molecular weight and strength of polar interactions, and the 
macroscopic rheology. Computer simulations can help in increasing this understanding, but 
no particle-based models are available that can reach large enough time and length scales 
to predict the experimentally accessible rheology. 

In this paper we introduce a novel highly coarse-grained simulation model that can predict 
semi-quantitatively the bulk rheology of a PSA. We will show that inclusion of transient 
forces {21, i.e. forces with memory of previous configurations, is essential to correctly predict 
the linear and nonlinear properties. Without transient forces even the qualitative features 
are incorrect. 

A large fraction of industrial waterborne is based on acrylic polymers synthesised 

via radical emulsion polymerisation processes [3]. Usually the formulation is complex, em- 
ploying random copolymers of a long-chain acrylic (such as n-butyl acrylate) characterised 
by a low glass transition temperature (Tp) with a short side-chain acrylic (such as methyl 
acrylate) to adjust Tg and acrylic acid to improve adhesion and optimise elongational prop- 
erties 0, S]. When the emulsion is dried, the latex particles coalesce and an adhesive film is 
formed. 

At first sight the adhesive film appears to be homogenous. However there are indications 
that a fundamental heterogeneity is retained within the structure, which is linked to an in- 
complete coalescence of the latex particles 6|, |7(|, see Fig. [TJ Memory of the original interface 
between latex particles could be retained for a number of reasons: (i) because surfactants, 
used during the emulsion polymerisation process and used to stabilise the latex dispersion, 
do not drain away fully but partly get stuck in the interfacial regions Q , (ii) because poly- 
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FIG. 1. (color online) Cartoon of the creation of an adhesive film. Left: polymeric latex par- 
ticles (white) are dispersed in water (blue) by the action of surfactants and polar groups which 
predominantly reside on and in the interfacial region of the latex particles (yellow). Right: when 
the system is dried, the water evaporates and most, but probably not all, surfactants are drained 
away. The polymer chains from different latex particles start to intermix, and polar groups cause 
a strong bonding between the particles. 



mer chains on different latex particles do not have enough time to fully interdiffuse, or (iii) 
because polar groups (such as groups of acrylic acid and sodium ions) preferably reside in 
the outer regions of the original latex particles and do not homogenise after coalescence. A 
combination of these effects is also possible. Heterogeneities in the concentration of remain- 
ing surfactants can be visualised by transmission electron microscopy, and heterogeneities 
in polymer and polar group concentrations can be characterised by atomic force microscopy 



|9|,|l0|. The effects of heterogeneities can also be traced back in the dynamical behaviour 
of the film. The bulk rheology for example is more complex than that of a homogeneous 
melt of polymers with comparable molecular weight and branching frequency. 

Not surprisingly, the design of new waterborne PSAs is a complex formulation process, in 
which the number of parameters is so large that the approach is often empirical and far from 
optimised. To enable a knowledge-based design, a fundamental understanding is needed of 
the link between the molecular properties and structure of the latex particles on the one 
side, and the mechanical properties of the adhesive film on the other side. Particle based 
computer simulations could provide such a link. 



Here we present a novel mesoscopic Brownian Dynamics (BD) simulation model in which 
the smallest unit is one full latex particle produced by emulsion polymerisation. This degree 
of coarse graining is as large as possible to permit a large integration step (3 10^^ s) and few 
degrees of freedom. The key issue is to identify the most important microscopic properties of 
a PSA that determine its rheology. In our work we have assumed that these key properties 
are the mechanical properties that characterise a single latex particle, such as its size and 
deformability, and the properties of the interfacial region, in particular the interaction free 
energies and relaxation times associated with latex interparticle adhesion and polymer chain 
intermixing. It is important to note that in our one-particle approach the latex particles keep 
their individuality, even after drying. It is therefore implicitly assumed that no major mass 
transfer between the latex particles will occur. Whether or not this is the case for real latex 
particles is presently unclear, but the observed heterogeneities suggest that interdiffusion is 
probably not dominant. 

Before embarking on a programme in which the model is tweaked and optimised to 
quantitatively predict the rheology of specific PSAs, our first goal is to test whether our 
model can predict the quahtative features. We found that using crude estimates for the 
model parameters already enabled us to reproduce semi-quantitatively the linear rheology 
of a particular experimental PSA, which we refer to as our "reference material". We will 
subsequently compare predictions for the non-linear extensional rheology with experimental 
measurements. As already mentioned, inclusion of memory effects, through application of 
transient forces, is necessary to correctly predict the linear and nonlinear dynamic properties. 

This paper is organised as follows. The BD model, including the treatment of transient 
forces, is explained in section [III Details of the reference material, its preparation and, cou- 
pled to this, a justification for the simulation parameters, are given in section UTTl In section 
IIVI we describe the rheological experiments we have performed on the reference material. In 
sections |V] and |Vl] we present the simulation predictions for the linear and nonlinear rhe- 
ology, respectively, and compare with experimental measurements on the reference system. 
We end with conclusions and an outlook on further developments of the model in section 

Em 
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II. MESOSCALE MODEL WITH TRANSIENT FORCES 



A. Coarse graining 

Our goal is to describe the configurations and forces in a pressure sensitive adliesive 
formed by drying of a latex emulsion by as few variables as possible, while at the same time 
retaining a link with the chemistry and individuality of the latex particles. We choose to 
represent the configuration by just six coordinates for each latex particle, corresponding to 
a vector locating its centre of mass position and a vector characterising its extension. This 
is not to say that all the removed coordinates are irrelevant for the rheology of the system. 
They provide an effective potential, a so-called potential of mean force A"^-^, governing the 
equilibrium distribution of the centres of mass and extensions. In thermodynamic equilib- 
rium the probability distribution P^q [r^^ , e^^^ of the particles' centre-of-mass positions and 
extension vectors is given by 

^m/ ^^3Af^ g37V^ 



Pe. p^,e=^^]ocexp 



(1) 



It is possible to derive equations of motion for the retained coordinates, not only leading 
to a mean force (which is minus the gradient of the potential of mean force), but also to 



frictions and random forces 
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In most coarse-grained models of soft matter systems, 
these frictions and random forces have memory of the configurations the system has gone 
through in the recent and sometimes even the distant past [2]. In such cases a simple 
Brownian dynamics propagator with realistic mean forces and uncorrelated, fully random 
displacements without memory will not reproduce correct sequences of configurations of the 
retained coordinates. To circumvent the introduction of memory effects in friction forces 
and stochastic displacements we employ the Responsive Particle Dynamics (RaPiD) method 



18|. 



The idea behind the RaPiD method is to introduce a relatively small set of additional 
dynamic variables which keep track of the thermodynamic state of the eliminated coordi- 
nates for the given values of the retained coordinates. The deviations of these additional 
variables from their equilibrium values at the given configuration give rise to strong forces in 
addition to the thermodynamic forces deriving from the potential of mean force. While the 
additional non-equihbrium forces gradually fade away because the variables relax towards 



their equilibrium values, with every new time step new such forces are being generated be- 
cause the displaced coordinates dictate new equilibrium values for the variables. For this 
reason the additional forces are called transient forces. 

In the following subsections we will introduce our model for the potential of mean force 
and the transient forces. We will give the equations of motion by which we update the 
retained coordinates and additional variables characterising deviations from equilibrium. 

B. Potential of mean force 

The latex particles are relatively soft and deformable. In the simulation we allow the latex 
particles to deform and orient by representing each particle, labeled i, as a hemispherical 
cylinder with an internal finitely extensible nonlinear elastic (fene) spring of length 2ei, as 
shown in Fig. |2l The diameter cTj of the cylinder is defined as o"o at zero extension (when 
the particle is actually a sphere). With increasing extension the diameter decreases such 
that the volume of the hemispherical cylinder remains equal to IvtcTq. An analytical function 
for the diameter may be derived by solving a cubic equation using Cardano's rules, but in 
practice in the range of extensions sampled in this work (0 < 2ei < 2.5(7o) a good fit is given 
by the much simpler function 



with Oi = 1.59 and 02 = —0.61. 

Figure [2] also shows how we let the interaction between each pair of particles be a function 
of the closest distance dij between the line segments running through the cylinder axes: 



Here Tj and ej are the location of the centre and (half) the extension of dumbbell i, respec- 
tively. Note that the values of Sj and Sj are constrained to the interval [—1, 1]. In practice, 
we first solve analytically the unconstrained minimisation of the bivariate quadratic function 
dfj{si, Sj). If the resulting values of Sj and Sj are not inside the square represented by Si = ±1 
and Sj = ±1, then each of the sides of this square are checked. The latter corresponds to 
a situation in which an end-point of the segment is in closest contact with the other line 
segment. 




(2) 



din 



"I3 ~ 




(3) 



7 



FIG. 2. Each latex particle is represented by a hemispherical cylinder located at r, with an 
orientation and extension given by e (cylinder length 2e). A finitely extensible nonlinear elastic 
spring is associated with each extension. All interactions between pairs of particles i and j are 
functions of the closest distance dij between the line segments running through the cylinder axes. 
The diameter cjj of a cylinder decreases with increasing extension Cj in such a way that the volume 
of the hemispherical cylinder is conserved. 

We let the particles interact with a repulsive interaction y^''^^ and an attractive interaction 
ip'^, both of which scale with the average particle diameter aij = + (Xj). Physically, the 
repulsive interaction yj''^^ = ip'''^'P{dij/(Tij) results from the entropic cost of chain deformation, 
as well as more direct excluded volume interactions between monomers. The former sort of 
repulsion is relatively soft and important at intermediate distances between latex particles, 
whereas the latter is stronger and more important at close distances. The latex particles do 
not only repel each other, but also attract because the polar groups generate strong bonds 
between individual chains. In the model we assume that the free energy associated with the 
cohesion between a pair of latex particles changes from — e'^ to zero with increasing scaled 
interparticle distance through a smooth 'sticker' potential function ip^ = ip^{dij/aij). 

In detail, the total potential of mean force has the following form: 



(4) 




(5) 



(6) 
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FIG. 3. (color online) Potential of mean force, with the parameters used in this work. Main figure: 
repulsive (dashed red line) and attractive sticker (solid red line) pair interaction as a function of 
distance dy for aij = ctq = 400 nm. The total pair interaction (solid black line) has an attractive 
well around do. Inset: potential of the finitely extensible nonlinear elastic spring as a function of 
the total intraparticle extension 2ej. The maximum extension in this work is Rq = 1000 nm. 



cp%d,a) = -e'n^,{d,a) (7) 

(8) 



1 -tanh ( — ( - - 1 ) - 1 



The fene potential in Eq. ([5]) behaves as a linear spring with spring constant k at small 
extensions, while it diverges asymptotically at the maximum extention Rq. The repulsive 
potential is characterised by its softness parameter n and equal to a^''^'^ at d = a. For low 
values of n the repulsive interaction is, as desired, relatively soft in the range a < d < 2^^"a 
and quickly becomes stronger for d < a, where it diverges as r~^". The particles can 
therefore interpenetrate, but not so far that they can fully 'cross each other'. The repulsive, 
sticker and fene interactions are shown in Fig. [31 

Note that we have chosen, rather arbitrarily, a tangenthyperbolic function for the sticker 
potential of mean force. The corresponding function nQ{d,a) may be thought of as the 
average fraction of sticker groups still active between a pair of particles at distance d. The 
only requirement for this function is that it is smooth and decays from one to zero over 



a finite width w^. We have chosen its inflection point at do + w'^ (scaled by cr/ao), which 
ensures that almost all sticker groups are active when d < a. 

We will give estimates for the parameters appearing in the potential of mean force in 
subsection IIII B[ In principle the potential of mean force can also be obtained directly by 
means of an atomic force microscope by measuring the average forces between two latex 
particles held at various distances, although in practice this has not been done yet. 



C. Transient forces 



For the current system we identify two major sources of transient forces, namely devia- 
tions in the number of sticker groups and deviations in the amount of 'intermixing' of the 
chains. The latter is a topological effect: when the distance between two latex particles is 
suddenly changed, the relaxation from unfavorable polymer conflgurations to new favorable 
ones is slowed down by the connectivity and uncrossability of the polymer chains. We denote 
the fraction of intermixing of the chains belonging to particles i and j by the symbol nfj, 
where the 'e' stands for entanglement, with the understanding that any topological effect 

is included: these are not literally the entanglements treated in the Doi and Edwards tube 
□ 

model [19|, but rather a measure for intermixing between the chains belonging to different 
latex particles. Similarly, the symbol n^j is used for the actual fraction of sticker groups 
shared between latex particles i and j. We assume that the equilibrium number of sticker 
groups and the equilibrium amount of chain mixing between i and j are functions only of 
the distance between these two particles. 

It is possible to show 2| that, to lowest order, the transient forces may be thought of as 
originating from a penalty free energy function A^^"'""'^ which is quadratic in the deviations 
of the additional variables from their equilibrium values, i.e. 



A 



trans 



E 

Kj 
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^0 i^ij^ '-^i. 



n, 



(9) 



where, similar to the average number of sticker groups, the average amount of chain mixing 
between a pair of latex particles is taken to be a smooth general S-shaped function: 

1 



1 



tanh f — f — 



a 



(10) 



Here gives the characteristic width over which the chain mixing decays smoothly from 
one for d < a to zero for d > a. 
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The new variables a* and a*^ control the allowed equilibrium deviations in the number of 
sticker groups and chain mixing. Generally, larger values of a'^ and correspond to larger 
penalties for deviations, resulting in stronger transient forces. 



D. Equations of motion for the retained coordinates 

In the preceding subsections we have defined the potential of mean force and the 
penalty free energy A^^"-'^'^ associated with transient forces. The actual simulation proceeds 
by integrating (by a simple first order Euler scheme) the following stochastic equations of 
motion: 

Here A = A"^^ + ^4*''""'* is the total free energy of the system and and ^'^ are frictions 
on displacements of the centre of mass and internal deformations, respectively. We need a 
model to describe these friction functions. We will assume that the frictions depend linearly 
on the amount of chain mixing with neighbouring particles. In particular, we use: 



^l = a=^eJ:^J^^,^0id^,)^ (13) 

were C,e is the maximum friction generated by the mixed interface between a pair of parti- 
cles. For simplicity, we have made the frictions on the positions and the extensions equal. 
Although this can be relaxed in future work, it is expected that the two frictions will not be 
very different because a change in position generates roughly the same amount of relative 



motion as a similar change in extension e,. In Eq. f|T3|) we use a form J2j \l^ij^1){d)i rather 
than J^j'iT'ij, because this allows the friction contributions to go to zero smoothly at the 
cut-off distance. In practice we use large values of a'^ and therefore n^j will always be close 
to Uq, and the difference will therefore be small. 

In principle, to ensure a correct equilibrium probability density, also a 'drift' term 
d/dri{kBT 1^1) should be added to Eq. ffTTl) . and a similar term to Eq. f|T2|) . but in our 
simulations the densities are always so high that the gradients in friction are negligible 
compared with gradients in the free energy. 
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Finally, the random forces P(t) and i^{t) are related to the frictions through the 
fluctuation-dissipation theorem, resulting in: 

{r{t)r{t')) = ^is{t-t') (14) 

Si 

rimt')) = ^-^mt - 1') as) 

Here I is the unit tensor and 6{t) is the Dirac delta function. 



E. Equations of motion for the sticker and entanglement numbers 



Besides the equations of motion for the retained coordinates, we also need equations of 
motion for the sticker and chain mixing fractions, nf,- and n%. These numbers tend to relax 
towards the equilibrium values prevailing at the given particle configuration. It is important 
to realise that such a relaxation is not instantaneous, but requires a certain amount of time. 
This is exactly what gives rise to memory and transient forces. 

We update the sticker and intermixing fractions according to: 



d?4 

dt 
dnf- 



n. 



1 

2kBT 
2kBT 



no (dij) 

6{t - 1') 
6{t - 1') 



+ 9\t) 



(16) 
(17) 
(18) 
(19) 



Here and are the characteristic times with which deviations in the fraction of sticker 
groups and chain mixing, respectively, relax towards their equilibrium values. For simplicity, 
in this work the two characteristic times are (different) constants. It is likely that the 
relaxation times actually depend on the distance between a pair of particles (i.e. their 
overlap). This may be included in future work when the model is refined further. 

The amplitudes of the random terms g'^ and in Eqs. f|T8l) and f|T9|) are chosen such that 
in equilibrium the expected fluctuations in sticker and mixing fractions are given by: 



I]) 



a" 
ksT 



(20) 
(21) 
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These latter two equations give us a convenient way to check (in equihbrium) whether the 
equations of motion are integrated with a sufficiently small time step [l3|: if the measured 
fluctuations in sticker and mixing fractions do not agree with Eqs. f l20|) and f l2T|) . a smaller 
time step must be used. 



F. Simulations without transient forces 



To assess the effect of transient forces, we also include simulations without transient 
forces. All conservative interactions, i.e. all terms in the potential of mean force A"^-^ , are 
retained. The transient penalty free energy A'^^"'^'^, however, is set to zero. For the friction 
we use exactly the same expression as before, Eq. (fT3l) . but manually set nf^ equal to the 
equilibrium number nQ{dij). 



G. Boundary conditions and nonlinear flow 

As with other particle-based methods, our simulations are typically limited to systems 
containing several thousands of particles (in this work 10,000). To mimic bulk behaviour, 
we apply periodic boundary conditions, where a primary box is surrounded by copies of 



itself 



20l ]. When a particle leaves the primary box, its image will simultaneously enter the 



box through the opposite face. Therefore the shape of the box must be space filling. We 
choose a rectangular parallellipiped, with box sizes along x and y twice as large as the box 
size along z, = Ly = 2Lz, for reasons that we will discuss next. 

In this paper we apply two types of deformation, namely planar (oscillatory) shear flow 
and planar extensional flow. 

Planar shear flow is characterised by a flow field 



^0 7 0^ 



V r 








r. 



(22) 



Such a flow field can be simulated by applying Lees-Edwards periodic boundary conditions 



211 ]. and augmenting the equations of motion with a flow given by Eq. fl22l) . The motion 



of a box image is such that its origin moves with a velocity in the x-direction which is 
proportional to the ?/-coordinate of that particular box's origin. If a particle moves through 
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the lower y-face of the primary box, the replacing image particle at the upper y-face will 
not have the same x-coordinate, but one displaced by an amount 5x = 'jLyt (where 6x may 
be taken modulus L^). 

Planar extensional flow is characterised by a flow fleld 



V r 



e 
-e 
\0 / 



(23) 



It is more diflicult to apply a continuous extensional flow to a periodic system. We use a 



method proposed by Todd and Daivis 
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231 ] and by Baranyai and Cummings 2J]. Both 



methods are based on work of Kraynik and Reinelt 25|. The basic idea is to choose a 
primary box which is square in the xy-plane and rotated over some angle 6. This angle 
is not arbitrary, but must be chosen such that after a flnite simulation time the deformed 
box coincides with some unit cell of the lattice generated by the original box. At this time 
all particles are mapped back into the original unit cell and the run is continued. Out of 
the inflnite set of solutions, we choose the case for which the minimum distance between 
opposite faces in the deformed box is still as large as possible. Nevertheless, for this case 
the minimum distance is a/5 times smaller than the original box size (or Ly) To 
prevent particles from interacting with their own images it is important that the box size is 



20|. This 



always at least twice the maximum interaction range between a pair of particles 
is the reason why the box dimensions in the x and y direction were chosen twice as large as 
in the z direction. 

The response of the system to flow is measured by keeping track of the instantaneous 
microscopic stress tensor for a periodic system 



Tj) Fij. 



(24) 



Here V is the volume of the primary simulation box, rj — is the vector pointing from 
the centre-of-mass of particle j to the centre-of-mass of particle i, and Fij is the total 
interparticle force (both conservative and transient) exerted on particle i by particle j. 
Because of the dominance of these interparticle interactions, we are justifled in neglecting 
any kinetic contributions to the microscopic stress. This is also implicit in our Brownian 
Dynamics approach. 
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III. THE REFERENCE MATERIAL 



In this paper we aim to model the rheology of a particular pressure sensitive adhesive, 
the reference material. In subsection IIII Al we will briefly introduce the reference material 
and its preparation. In subsection IIII Bl we give an overview of the parameters used in our 
simulations and a justification for their values. 

A. Materials and preparation 

The reference material consists mainly of acrylic polymers with a mean molecular weight 
of 2203 kg/mol. A latex emulsion with 53.6 % solid content was received from DOW Chemi- 
cal. An adhesive film was prepared by uniformly depositing about 2 ml of adhesive emulsion 
into a fiat rectangular silicone mold cavity (23 x 47 x 6 mm), which was then allowed to air 
dry at room temperature for about one week. Finally, the now transparent films were dried 
in a vacuum oven for 5 minutes at 80°C to remove remaining water. These dried films were 
removed from the mold cavity and sandwiched between silicone papers in order to avoid 
sticking to any other substrate during handling. 

B. Choice of simulation parameters for the reference material 

An overview of the simulation parameters and the particular values used in this work is 
given in Table HI Most parameters of the model have been estimated from available data 
(T, 0, (To, w'^ and r"^) or statistical reasoning (n, a*"*^^, k, Rq, e**, w"^ and ^'^). A few parameters, 
(r*, a** and a^) could not be estimated and were treated as fit parameters to reproduce the 
experimental linear rheology. 

Transmission electron microscopy has revealed that before coalescence the latex particles 
have a diameter of 400 nm. After coalescence and drying of the film most solvent has 
evaporated. Assuming that the particles have not significantly swollen or shrunk during 
the drying process, we have set our bare particle diameter to ctq = 400 nm. We have 
chosen the number density p of latex particles such that the corresponding volume fraction 
(j) = Ipcj^ is equal to 1.0. For hard spheres such a high volume fraction would be unachievable 
{(pmax ~ 0.64 for random close packing) but it should be remembered that the latex particles 
have a polymeric content and are therefore soft and deformable. This softness is reflected 
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property 


symbol 


value 


General parameters 


temperature 


T 


303 K 


voluuw^ fraci ion 


6 


1.0 


Particle 


diameter 




400 nm 


softness parameter 


n 


6 


repulsive energy (pair) 




500 ksT 


internal spring constant 


k 


3 10-3 N/m 


maximum extension 


Rq 


1000 nm 


Sticker clusters 


adhesive energy (pair) 




4 10^ ksT 


transition width 




40 nm 


transient penalty 


a' 


8 10^ ksT 


relaxation time 




1 s 


Interparticle chain mixing 


transition width 




40 nm 


transitient penalty 




10^ ksT 


relaxation time 




10-3 s 


friction per particle pair 




2 10-"^ kg/s 



TABLE I. Simulation parameters representing a reference pressure sensitive adhesive based on 
acrylic polymers. The interparticle chain mixing time and transient penalties a* and are fit 
parameters to semi-quantitatively tune the linear storage and loss modulus. All other parameters 
have been estimated from available data or statistical reasoning. 

in a low value, n = 6, for the scaling of the repulsive potential. We have also chosen a 
relatively low value for the repulsive energy parameter a^^^ which is discussed after treating 
the sticker forces. Together with the attractive sticker forces, the particles can interpenetrate 
quite easily (in contrast to hard spheres), with closest distances as small as d — O.Tctq and a 
peak in the radial distribution function at d = O.SSctq (not shown). 

The cohesion between the latex particles is mainly caused by polar charge interactions in 
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clusters of acrylic acid -COOH or sulphate -OSO3 centered around sodium or other positively 
charged metal ions. The binding energy for two carboxylates with identical sodium (in the 



absence of water) is of the order of 40 fc^T 26|, |27|. Each latex particle contains a very large 



number of polar charge groups on its surface. Based on the preparation chemistry, and the 
fact that each latex particle on average interacts with 10 neighbouring latex particles, we 
estimate that there are typically 1000 polar bonds between each pair of latex particles. This 
leads us to an estimate of e'^ = 4 10^ ksT for the total adhesive energy between a pair of 
latex particles. 

The characteristic width over which the attractive "sticker" potential <y9'^ changes from 
— e'^ to zero is not the distance over which a single polar bond breaks. It is rather the distance 
over which the number of polar bonds shared between two latex particles changes from its 
maximum to zero. This will depend on the precise distribution of polar groups through 
the interfacial region of the latex particles. Unfortunately this information is lacking. We 
therefore assumed that the polar groups are scattered more-or-less randomly throughout the 
interfacial region of the latex particle. Changes in the number of shared polar groups should 
then occur over the same length scale as changes in chain intermixing, i.e. = w^. The 
chain intermixing will typically change on a length scale equal to the radius of gyration of 
the polymer backbone. For a closely related = 2203 kDa PMMA chain, the radius of 



gyration is estimated to be about 37 nm 28|]. In the model we use w'^ = w'^ = 40 nm, which 
is 10% of the particle bare diameter (Tq. 

Given the above estimates for and w^, we choose the repulsive energy parameter a^'^^ 
such that the attractive well (f^^^^ + (f^ is more or less symmetric, with a characteristic width 
of 2w^, as shown in Fig. [31 

At this point the potential of mean force of our model is almost fully specified, except for 
the parameters of the finitely extensible spring, which determines the internal deformability 
of a latex particle. 

The range of maximum extension is set to Rq = 1000 nm. Changing the value of Rq 
hardly has any influence on the linear rheology, but it can change the nonlinear rheology. 
For the purpose of this paper, it is important that Rq is considerably larger than the particle 
bare diameter ao so that the particle can accommodate a large stretch under the influence 
of strong nonlinear flow. The exact influence of Rq on the nonlinear rheology is a topic of a 
future paper. 
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Next we must choose the spring constant k. When our model material is under a stress, 
it can relax this stress by changing interparticle distances and/or particle extensions. If the 
force associated with changing a particle extension over a certain distance is much weaker 
than the force associated with changing the interparticle distance over the same distance, the 
system will always respond by elongating the particles. In the opposite case, the system will 
always respond by changing interparticle distances, while keeping the particles effectively 
spherical (elongation zero). In real PSA materials these two effects are more or less balanced, 
so it is important to choose the value of the spring constant k roughly equal to the curvature 
of the total interparticle pair interaction near its minimum. This way we arrive at A; = 3 10~^ 
N/m. 

The friction per particle pair was chosen large enough to suppress Brownian motion 
(which is also suppressed in the real material) and allow for a large integration time step of 
3 10~^ s, but not so large that interparticle distances and extensions are effectively frozen 
on the time scales of interest (10^'^ to 10 s). A good compromise was found at = 2 10~^ 
kg/s. 

Now we need to specify further the parameters for the transient forces. Consider again the 
case where the distance between a pair of neighbouring latex particles is suddenly increased. 
Even if the polar groups would be very weak and break easily, the intermixed chains in the 
interfacial region between the two latex particles would need a finite time to relax to new 
configurations. Typically, such a relaxation takes place on the Rouse tinae scale, which for 
our acrylic chain has been measured to be of the order of a millisecond [29']. We therefore 
choose = 10~^ s. Of course the polar groups also need a finite time to break, and in view 
of the large binding energy associated with a polar bond, it comes as no surprise that this 
time should be relatively large. Unfortunately no measurement or estimates are available, so 
currently we treat this as a fit parameter. For our reference material we find good agreement 
for the dependence of the loss modulus on frequency if we choose r'* = 1 s, which seems 
reasonable. 

The transient penalties a** and have been used to semi- quantitatively tune the simu- 
lated storage and loss modulus to the experimentally measured ones in the range u = 10~^ 
rad/s to 10^ rad/s. We found that, given the relative smallness of t'^, the parameter a'^ only 
influenced the linear rheology at w > 10^ rad/s. Therefore, a* was varied first to tune the 
linear rheology for low frequencies u < 10^ rad/s, leading to a"^ = lO^ksT, after which 
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was varied to also tune the high frequencies, leading to — lO^/csT. 

In a forthcoming paper we will further investigate the influence of each of the model's 
parameters on the rheology. Here we already note that the parameters which are expected to 

be most strongly dependent on the PSA chosen for study are e'*, ^ and together char- 
acterizing the distribution and strength of the sticker groups, and w^, and i?o, together 
characterizing chain entanglement and stretch. 

IV. EXPERIMENTAL METHODS 

Disks of 8 mm diameter were cut out of the dried films by a precise punching tool for 
the shear flow tests, while for the extensional flow characterization, rectangular sized films 
of 1-2 mm width and about 15 mm length were prepared using a razor blade. 

For the linear viscoelastic characterization, small-amplitude oscillatory shear-flow tests 
were performed with a strain-controlled rheometer (ARES, TA Instruments) using parallel- 
plate tools with diameters of 8 mm at frequencies from 10~^ rad/s to 10^ rad/s and at 
temperatures from -45°C to 90°C. These results were used (through time-temperature su- 
perposition) to extend the range of accessible frequencies at the reference temperature of 
30°C. 

The nonlinear- viscoelastic behavior of the PSA was measured in uniaxial elongational flow 
using the extensional viscosity flxture (EVP, TA instruments) on the same rheometer at a 
temperature of 30°C and at strain rates of 1 s~^ and 10 s~^. Elongational stress relaxation 
was measured by cessation of extensional flow at strain rates of 1 s~^ for different Hencky 
strains. 

V. LINEAR RHEOLOGY 

We start by exploring up to what point our model PSA behaves hnearly. The material 
is subjected to an oscillatory shear, with the flow in the x-direction, a gradient in the y- 
direction, and a strain given by 

7(t) =7osin(a;t), (25) 

where 70 is the strain amplitude and uj the oscillation frequency. In the linear regime, 
the shear stress response is related to the storage modulus G'iiS) and loss modulus G"{yi) 
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FIG. 4. (color online) Strain amplitude sweep test of the linear regime for the modeled PSA under 
oscillatory shear at w = 1 rad/s. Shown are the storage modulus G' and loss modulus G" (closed 
and open squares, left scale) and the loss tangent tan 6 (triangles, right scale) The material behaves 
linearly (constant G' , G" and tan 5, dashed lines) up to a strain of approximately 0.1. Experimental 
results, included for reference (green symbols), are linear up to a strain of approximately 0.4. 

according to 

Sxyif) = loG'i^) sin (out) + '^qG"{uj) cos {ujt) . (26) 

Outside the linear regime, the stress response is no longer linear in 70 and possibly higher 
harmonics (Sw, Sw, . . .) will appear. In all cases, the base frequency responses G'{uj) and 
G"{uj) have been calculated from 

G'{uj) _ J_ /■* ^ (t) sin(wt)dt, (27) 
7or Jt' 

G"{uj) _ J_ /■* 5 (t) cos(a;t)dt, (28) 
7or Jt' 

where r is an integer number of oscillation periods (r = 27171 /u, with n an integer), and the 
starting time t* was chosen sufficiently large to avoid any startup effects. 

For strain amplitudes up to approximately 0.1 our model behaves linearly, as evidenced 
by the constancy of G' and G" as a function of 70 in Fig. H] (black symbols) and the absence 
of higher harmonics in the stress response (not shown). Beyond this value higher harmonics 
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FIG. 5. (color online) Linear storage (filled symbols) and loss (open symbols) modulus. The main 
figure shows the results of simulations with (black squares) and without (red squares) transient 
forces, as well as experimental results of the reference system (green circles). The inset shows where 
our targeted range of frequencies fits in the overall range of experimental frequencies, as obtained 
by time-temperature superposition for a reference temperature of 303 K. 

emerge and both the storage and loss modulus start to drop, but the storage modulus drops 
faster. As a consequence the loss tangent tan 5 = G"/G' increases slightly beyond a strain 
amplitude 0.1. Experimentally, the linear regime extends up to a strain amplitude of 0.4 
after which similar nonlinear deviations are observed, as shown by the green symbols in 
Fig.H 

In our simulations, the linear storage and loss modulus have been determined for frequen- 
cies ranging from 10~^ rad/s to 10^ rad/s. In this range of frequencies the storage modulus 
G' is consistently larger than the loss modulus G", hence the material's response is mainly 
elastic, see the main plot of Fig. |5l This is the case both with and without transient forces. 
Without transient forces, however, the dependence of the loss modulus on frequency is rather 
strong and Maxwell-like (open red squares in Fig. [5]). The influence of the transient forces, 
with the current parameters settings, is to increase the (elastic) storage modulus at higher 
frequencies and, more significantly, to increase the (dissipative) loss modulus at both high 
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and low frequencies (open black squares). As a consequence both the storage and loss mod- 
ulus increase monotonically with increasing frequency. The increase is slow, spanning less 
than a decade compared to four decades increase in frequency. The increase is particularly 
slow for the loss modulus G" for frequencies up to 10 rad/s. 

Fig. 13 (main plot) shows that qualitatively all the above features in the storage and loss 
modulus are also present in the experimental system. Given the crude way in which we have 
constructed our mesoscale model, the near-quantitative agreement with the experimentally 
determined storage and loss modulus is encouraging. We expect that only moderate re- 
finements to the model and its parameters will be needed to achieve a full quantitative 
agreement. Such refinements will be different for each particular pressure sensitive adhesive; 
the dependence of the rheology on the simulation parameters will be explored in a future 
paper. 

For reference, in the inset of Fig. [5] we show the linear rheology of the system in a much 
wider range of frequencies, 10~^ rad/s <uj < 10'' rad/s, as can be obtained experimentally 
by time-temperature superposition. Because our particles are coarse-grained, we do not 
expect nor do we aim for agreement at frequencies higher than 10^ rad/s. On the other side 
of the spectrum, because of computational limitations, our simulations are currently limited 
to a total run time of order 10 seconds. We can therefore not yet access the w < 0.1 rad/s 
regime. With further increase in speed of computer hardware, combined with parallellization 
and optimization of the code, we expect to be able to reach one or two orders of magnitude 
lower frequencies in the near future. 

VI. NONLINEAR EXTENSIONAL RHEOLOGY 

We will now test the predictions of the mesoscale model in nonlinear flow. One of the 
most extreme forms of nonlinear flow is extensional flow. Unlike shear flow, extensional 
flow has no rotational component, leading to a continual relentless extension of the particles 
without a possibility to escape through reorientation. 

As described in section III G\ we apply a planar extensional flow with constant extension 
rate e to a previously equilibrated sample. This flow is suddenly started at t = and 
the time-dependent extensional stress Ni{t) = S^xit) — Syy{t) is measured. From this we 
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FIG. 6. (color online) Extensional viscosity growth rj'^{t) as a function of time t after a sudden onset 
of a planar extensional flow, for two different extension rates e (see legend). We show results for 
simulated systems with (black lines) and without (red lines) transient forces. Extensional viscosity 
growth predicted from the linear rheology (see main text for an explanation) is shown as a blue line. 
In the experiments a uniaxial extensional flow is applied; the results for the extensional viscosity 
growth are shown for reference (green symbols). 

calculate the time-dependent extensional viscosity, defined as 



In Fig. in] we show our results for two extension rates, e = 10 s~^ (black solid line) and 1 
s~^ (black dashed line). In the following we will compare these results with (i) a prediction 
based on the linear rheology, (ii) results from simulations without transient forces, and (iii) 
experiments. 

For Newtonian liquids the extensional viscosity is related to the shear viscosity through 
a multiplicative constant, referred to as the Trouton ratio (Tr). For planar extensional flow 
Tr= 4, while for uniaxial extensional flow Tr= 3. Similarly, for non-Newtonian liquids the 
time-dependent extensional viscosity may be compared with the following linear prediction: 



(29) 



e 




(30) 
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where 'r]*{u}) is the complex shear viscosity as obtained from hnear oscillatory shear mea- 
surements. This prediction, shown as a blue line in Fig. El is in good agreement with the 
measured extensional viscosity at small times. At larger times, however, nonlinear effects 
start to dominate and the extensional viscosity grows faster than predicted from the linear 
model. For both extension rates the strain hardening starts at a Hencky strain en = of 
approximately 0.3. 

There are clear qualitative differences when comparing simulations with and without 
transient forces. In going from lower to higher extension rate, the order of the magnitudes 
of the viscosity are reversed: without transient forces (red curves), higher extension rates 
result in lower extensional viscosities at any particular time beyond the linear regime. Thus 
without transient forces, the simulations predict strain thinning behaviour. We can conclude 
that for this material, the occurrence of strain hardening behaviour depends crucially on the 
inclusion of transient forces. 

Currently no experimental data is available on planar extensional flow of the PSA ref- 
erence material. However, we can make a qualitative comparison with the extensional vis- 
cosity measured experimentally in uniaxial extensional flow. These results are shown as 
green symbols in Fig. El It is encouraging to observe that a similar strain hardening occurs 
experimentally, although the onset is somewhat later, at a Henky strain of approximately 1. 
We expect a similar onset of nonlinear effects in experimental planar extensional flow. The 
approximately threefold increase in the linear strain regime is consistent with the observed 
threefold increase in the linear regime under oscillatory shear described in section |Vl 

At a certain Henky strain a maximum extensional viscosity is reached. Experimentally, 
this occurs at a higher Henky strain than in our simulations. This could be related to the fact 
that the maximum extension of the latex particles in the simulations is lower than in reality, 
leading to an extensional stress which reaches its maximum sooner. We test this hypothesis 
in Fig. [71 which shows that indeed the mean square extension is close to the maximum 
extension at the time when the extensional viscosity reaches its maximum. The real 
polymeric particles can be extended to a larger amount, leading to larger stresses and larger 
maximum strain values. In future work we will vary the value of Rq and investigate its 
influence on the nonlinear rheology. 

Finally, we study the relaxation of extensional stress after a sudden stop of extensional 
flow. This is a sensitive probe of the relaxation processes that take place as the system 
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FIG. 7. (color online) Mean square extension ((2e)^) of the latex particles as a function of time t 
after a sudden onset of a planar extensional flow, for two different extension rates e (see legend). 
The results are normalized by the square of the unperturbed particle diameter, cjg. We show results 
for simulated systems with (black lines) and without (red lines) transient forces. The maximum 
extensibility used in this paper, Rq = 2.5cjo, is indicated by the horizontal dotted line. 

evolves from a highly non-equilibrium to a new equilibrium state. Figure [8] shows the stress 
when a sample is (planar) extended with an extension rate e = 1 s~^, which is suddenly 
stopped at a Hencky strain of 0.33 (black dashed line), as well as the stress when a sample is 
extended with an extension rate e = 10 s~^, which is suddenly stopped at a Hencky strain of 
1.0 (dashed black line). We observe that the extensional stress has a small initial relaxation, 
which is faster when the Hencky strain is higher. Soon hereafter we encounter an extremely 
slow extensional stress relaxation. 

Using our experimental equipment, it is difficult to measure extensional stress relaxation 
after a high initial extension rate of e = 10 s^^. It is however possible to measure it after 
extending with a rate of 1 s~^. Figure [8] shows these results for two experiments, in which 
uniaxial extensional flow is stopped after a Hencky strain of 1 (green circles) or 2 (brown 
circles), respectively. The experiments conflrm the simulation predictions of a fast initial 
relaxation which is faster for higher Hencky strain. Furthermore, the rate of the very slow 
stress relaxation that follows after the initial fast decay is very similar to what our simulations 
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FIG. 8. (color online) Relaxation of extensional stress after a sudden stop of planar extensional 
flow, for an initial extension rate of e = 1 s~^ stopped at t = 0.3 s (dashed line) and an initial 
extension rate of e = 10 s^^ stopped at t = 0.1 s (solid line). Results are shown as a function of 
Hencky strain en = as calculated from the initial extension rate. In the experiments a uniaxial 
extensional flow is applied; results are shown for the extensional stress relaxation following an 
initial extension rate of e = 1 s~^ which is stopped att = lsort = 2s (green and brown circles, 
respectively) . 

have predicted. 



VII. CONCLUSION AND OUTLOOK 



We have shown that the rheology of pressure sensitive adhesives can be reproduced nearly 
quantitatively by a coarse-grained model if, besides conservative forces, also transient forces 
between the constituent latex particles are taken into account. The conservative forces 
consist of attractive interactions caused by polar sticker groups, a soft repulsion and a 
nonlinear elastic spring associated with deformations of the latex particle. The transient 
forces arise from temporary deviations in the number of stickers and chain intermixing 
between pairs of latex particles. 

Our model reproduces nearly quantitatively the dominance of the storage over the loss 
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modulus, their slow and monotonous increase with frequency, and predicts nonlinear elon- 
gational effects such as the onset of extensional strain hardening and the characteristics of 
the relaxation of the extensional stress after cessation of extensional flow. 

Further improvements of the model and tweaks of the parameters are possible. For 
example, the relaxation times r** and r*^ are now taken as constants, but may actually 
depend on the distance between a pair of particles. This will effectively generate an even 
wider spectrum of relaxation times. Also the friction on the centre-of-mass coordinate of 
each particle is now taken to be the same as the friction on the extension coordinate of 
a particle, and both are based on the amount of intermixing (overlap) with neighbouring 
particles. These frictions may be (slightly) different, and may also depend on the number of 
stickers shared with neighbouring particles. Finally, the value of the maximum extension of 
a latex particle is important in determining the strain (both in shear and extension) at which 
nonhnear effects in the rheology start to play a role. We expect that such tweaks to the 
model and its parameters will bring the rheological predictions in quantitative agreement 
with experimental results. 

The purpose of this paper was to introduce the mesoscale model and to show the im- 
portance of including transient forces. The qualitative and near-quantitative agreement 
with experimental results motivates us (i) to further investigate the influence of each of the 
model's parameters on the rheology, (ii) to subject the model to more complex types of flow, 
and (iii) to study (onset of) flow-induced instabilities such as cavity formation. Work along 
these lines is currently in progress. 
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